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Dynamics of Q-Balls in an expanding universe 
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We analyse the evolution of light Q-balls in a cosmological background, and find a number of 
interesting features. For Q-balls formed with a size comparable to the Hubble radius, we demon- 
strate that there is no charge radiation, and that the Q-ball maintains a constant physical radius. 
Large expansion rates cause charge migration to the surface of the Q-ball, corresponding to a non- 
homogeneous internal rotation frequency. We argue that this is an important phenomenon as it 
leads to a large surface charge and possible fragmentation of the Q-ball. We also explore the devia- 
tion of the Q-ball profile function from the static case. By introducing a parameter e, which is the 
ratio of the Hubble parameter to the frequency of oscillation of the Q-ball field, and using solutions 
to an analytically approximated equation for the profile function, we determine the dependence of 
the new features on the expansion rate. This allows us to gain an understanding of when they 
should be considered and when they can be neglected, thereby placing restrictions on the existence 
of homogeneous Q-balls in expanding backgrounds. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Q-balls Q are non-topological solitons that exist as 
minimum energy configurations of a complex scalar field 
(typically a flat direction in a MSSM model 0, S 0) 
with a global U(l) symmetry; the gauged version has 
been studied in 5]. First introduced by Coleman 0, 
following pioneering work by Friedberg, Lee and Sirlin 
, their formation, particularly through an Affleck-Dine 
condensate collapse 0, H, @ , has been frequently inves- 
tigated . However, little work has gone into looking at 
the behaviour of general fully formed Q-balls in a cos- 
mological background. The need to do this is becoming 
increasingly important as Q-balls are often discussed in 
a cosmological context, typically as a possible catalyst 
for baryogenesis Q and as candidates for dark matter 
0, . In this short paper we begin this investigation by 
considering light Q-balls (where the Q-ball energy den- 
sity is dominated by the background energy density) in 
a cosmological setting. 

The background expansion leads to changes in the Q- 
ball behaviour from the static case. We identify a pa- 
rameter e, which is the ratio of the Q-ball oscillation 
frequency (in the core of the Q-ball) and the Hubble pa- 
rameter, as a measure of the deviation from static back- 
ground solutions. We then show that in an expanding 
background the Q-ball characteristics are determined by 
e. By numerically obtaining the explicit relation between 
e and the deviations from the static solutions, we are 
able to answer two important questions: when should we 
be considering background expansion effects and how do 
these effects depend on the rate of expansion? Turning 
explicitly to the effect of the expansion on the Q-ball we 
show that the Q-ball parameters exhibit expansion driven 
oscillations and that the expansion causes the field rota- 
tion frequency to become non-homogeneous, driving the 
conserved charge to the surface of the Q-ball but not ra- 
diating it, whilst maintaining a constant physical radius. 
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In what follows, we begin by summarising the basics of 
Q-balls in flat space leading to a discussion of the prob- 
lems a dynamic background introduces and our proposed 
way of analysing the new solutions. In section ILTT1 we nu- 
merically solve the full system of a light Q-ball in an 
expanding Universe, and explicitly demonstrate the con- 
sistency of introducing the parameter e. We introduce 
our analytic approach in section \N\ and discuss its uses 
and implications in section We summarise our results 
in section IvH 



II. Q-BALLS IN AN EXPANDING 
BACKGROUND 



The evolution of a complex scalar field ip with a La- 
grangian density £ = \dip\ — V(\p\ 2 ) is given by 
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a spherically symmetric field of the form 
ip = e f(r), where the frequency lu and profile function 
/ are real, in a flat background we have 
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where V e ff = \u? f 2 — \V(f 2 ). Equation (|2.2I) is equiva- 
lent to Newtonian motion, with a time dependent friction 
term, of a particle in a potential V e ff, with r acting as 
time. A Q-ball corresponds to a solution where the par- 
ticle starts at some point ipo = ip(r = 0) and rolls down 
the potential coming to rest at the origin. This would be 
a field configuration with p — > as r - » oo, a necessary 
condition for finite energy, if we require V(0) = 0. The 
conditions for such a solution to exist are that V e ff is 
decreasing near the origin (allowing the particle to slow 
down and come to rest there) and that at the starting 
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point V e ff ^ | cpo 1 2 ^) > 0, so the particle has enough energy 
to reach <p — 0. Hence for a flat space Q-ball solution to 
exist, the value of uj is restricted to be in the range 
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Assuming appropriate values of w are found, then equa- 
tion l|2.2l) has an exact soliton (Q-ball) solution with a 
stabilising (against decay) Noether charge given by pj 



Q=lu \ip\^d 6 x 
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In a Robertson- Walker universe with scale factor a(t), 
the equation of motion becomes 
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du 



(2.5) 



In this case the ansatz ip = e luJt f(r) with constant uj is 
no longer a solution and so the flat space Q-ball scenario 
of a spherically symmetric field configuration with a con- 
stant internal rotation frequency does not hold. Now, in 
most cosmological settings we may neglect the Hubble 
drag term as a first approximation, however a better and 
more instructive approach is to keep all the terms and 
simply work in conformal time rj defined by dt = a{rf)dr). 
Introducing (f> = ipa(rj) we obtain the following evolution 
equation 



(2.6) 



where time derivatives with respect to 77 are de- 
noted by a prime. The utility of this approach is 
that all the time dependence of the background has 
been re-written as a time dependent potential for <f> 1 
VM 2 ) = a*V(\<j)\ 2 /a 2 )-(a"/a)\<j>\ 2 . 



III. NUMERICAL SIMULATIONS 



For our analytic approach in section HVI we introduce 
a small parameter e = H/loq, so for our purposes a natu- 
ral background for looking at the evolution of the Q-ball 
profile is de-Sitter space. In this case the parameter e is 
a constant and so by looking at the profile evolution over 
a constant range of the scale-factor (constant number of 
efoldings), but at different rates of expansion, we can ex- 
plore the dependence of the approximation on e. First of 
all though we shall perform simulations of the full evolu- 
tion, (|2.ti|) . so we can later test our analytic method. In 
our simulations we used a form of the potential motivated 
from supersymmetry, namely 

V(tp) = cap 2 - p(p 4 + 7/ (3.7) 

with a = 2.3 x 10 7 GeV 2 , (3 = 1.0, and 7 = 1.0 x 
10 _8 GeV~ 2 . The chosen values correspond to the small 
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FIG. 1: Figure showing uj in the core of the Q-ball as a func- 
tion of the scalefactor. 



field expansion of a SUSY potential of the form V{<p) — 
M 4 ln(l + tp 2 /M 2 ) where M ~ 10 4 GEV is the scale of 
SUSY breaking. Recall that we are considering light Q- 
balls, defined as those whose dynamics are determined 
by the background spacetime, and whose own contribu- 
tion to the energy density in the Friedmann equation is 
negligible. Given that, we choose an appropriately large 
value H ~ 10 2 GeV for the Hubble parameter. This cor- 
responds to a background energy scale of 10 10 GeV, or 
an evolution in the radiation dominated era, just after 
say reheating at the end of a period of early Universe in- 
flation. In order to test the robustness of the parameter 
e in determining the true solution, we ran a series of sim- 
ulations in which we fixed the charge of the Q-ball and 
varied H , increasing it from a starting value of lOOGeV. 
We verified our results for a number of values of Q but 
the results we show here correspond to Q=300, (giving 
e = 0.05 for H = lOOGeV). 

Starting with the profile of a static configuration, we 
obtained the evolution of the Q-ball as we turned on 
the expansion of the universe, reaching a de-Sitter back- 
ground by a = 1. We then evolved the profile with this 
background and a fixed value of e for ln(25) efoldings. 
This was repeated for different values of e. Throughout 
the evolution we found that the Q-ball charge (defined by 
the Noether charge —Airi J dr r 2 [cj)*4>' — (</>')* <f\ where r 
is the comoving radial coordinate) was conserved. So the 
expansion did not lead to leakage of charge. 

Numerically the internal rotation frequency (in comov- 
ing units) was defined as 



oj(r) — Im 



(r) 



(3.8) 



and we observed that for slow expansion rates its value 
at the core of the Q-ball exhibited oscillatory behaviour 
about a general evolution where it remained proportional 
to a as shown in Fig. ^ Similarly, the profile amplitude 
<f> also uniformally increased as a whereas the Q-ball co- 
moving radius R decreased proportionally to a" 1 . Re- 
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FIG. 2: Figure showing the charge density radial dependence 
for a profile in an expanding background (e=0.25), with the 
inset showing a typical static background profile. Note how 
the expansion causes the charge to move towards the surface 
of the Q-ball, compared to the static case. 



calling the definition of <f> = (pa, we see that the physi- 
cal Q-ball parameters (radius and amplitude cp) oscillate 
about constant values. The main deviation from a flat 
space Q-ball arose through the non- homogeneity of u>(r), 
which emerged as e increased. Recalling that for homoge- 
neous rotation frequencies the charge distribution is pro- 
portional to the amplitude of the profile function, and 
that no major change to the profile shape was observed, 
such a non-homogeneous lo leads to a new charge density 
distribution. This is evident in Fig. [3 which shows how 
the expansion (with e = 0.25) causes charge to move to- 
wards the surface of the Q-ball as opposed to the case 
of a typical charge distribution in a static background 
indicated by the inset figure shown. However, once es- 
tablished, the charge density profile did not vary with the 
expansion. The oscillations of the Q-ball parameters are 
to be expected in an expanding background; they cor- 
respond to oscillations about the minimum energy con- 
figuration of the Q-ball and are driven by the expansion 
of the universe. To see this, recall that the expansion 
leads to time-varying potentials, in which the minimum 
is varying with time, whilst the field is oscillating about 
the minimum. 



IV. ANALYTIC APPROACH 

The usefulness of adopting conformal time emerges be- 
cause l|2.6[l looks like (|2.2[) but with a time varying poten- 
tial. It is natural, therefore, to make the approximation 
that a Q-ball profile at a particular time is given by solv- 
ing (|2.6|l at that time, (j> = e iujr, f(r), leading to 

^ + m(r ,,2+ £ ,2 -r' VU) ) = '>- (4 - 9) 



FIG. 3: Figure showing the numerical end profile (solid line) 
and the approximated profile (dotted line) obtained from the 
static solution for the case e = 0.1. 

This equation is of the same form as the flat space equa- 
tion (12.21) . but the price we pay is the introduction of a 
new time dependent effective potential. At a particular 
moment in time, equation l|4.9|l has soliton solutions with 
the same constraints as the flat space case, but these are 
Q-ball solutions in a dynamic background. This is an 
appealing result; we can think of Q-balls in a dynamic 
background as flat space Q-balls with a time varying po- 
tential. It is clear that our approximate ansatz is not an 
exact solution as the amplitude has no time dependence 
to match the evolving time dependent potential. Solving 
for static profiles at different values of a and a" provides 
an approximation for the evolving profile. Equation l|2.5l) 
approaches the flat space case in the limit — — > 0. There- 
fore we parameterise our approximation by the adiabatic 
parameter e = ^ where H is the Hubble parameter and 
u)q is the frequency in the core of the Q-ball at the ini- 
tial scalefactor value. We now show that we can closely 
predict the evolving profile of the true solution, based on 
the corresponding static solution. And that deviations 
from the static solutions are characterised by e. 

We would like to establish how in an expanding back- 
ground the deviations observed from the static profiles, 
depend on e. It is straightforward to solve for a static 
profile with the same overall charge as the exact solution 
and at the corresponding values of a and a" . Having done 
this we obtain an adiabatic approximation to the evolved 
profile, which can be compared with the true solution. As 
an example of the strength of the approximation for say 
e = 0.1, Fig. |21 shows an overlap of the true final profile 
(solid line) and the corresponding 'static' approximation 
(dotted line); they show good agreement. To quantify 
the accuracy of the approximation we introduce the pa- 
rameter 

y=i_ i/^r >0 (410) 

where (j>N{r) is the numerically evolved profile and 4>A{f) 
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FIG. 4: Figure showing lnXave as a function of -. 



is the adiabatic approximation. Note that \ measures 
how the profiles differ including the phases of the pro- 
files, so such measurements would include effects from 
non-homogeneous rotation frequencies. Calculating x 
throughout the evolution we find that it exhibits oscil- 
lations of the same nature as those found in the Q-ball 
parameters. Averaging over these oscillations gives us a 
parameter Xave &s a measure of how well the adiabatic 
approximation holds as a function of e. Figure0|shows a 
plot of In x&vc against - , with error bars corresponding to 
the standard deviation of the oscillations from the mean 
(note that these are not errors in the value of Xave but 
are there to show the scale of the oscillations). We see 
that there is a negative logarithmic correlation meaning 
that the adiabatic approximation rapidly improves with 
decreasing e. 

Figure 0] helps us to understand when we should con- 
sider expansion effects in Q-ball dynamics and how they 
vary with e in that it provides direct information on how 
the expansion affects the profile. For example, for ex- 
pansion rates of e = 0.2 we already have deviations in 
the profiles corresponding to around 10%, whereas for 
e = 0.05 these reduce to a negligible value of 0.1%. There 
is another useful non-trivial aspect to Figure 4. We, of 
course, expect that as e — > 0, we recover the static so- 
lution. However, an important question is how rapidly 
do we approach it? We can see that the deviation from 
the static Q-balls depend on high powers of e. Such a 
strong, highly non-linear, dependence is difficult is derive 
analytically and it provides important information in a 
cosmological setting where e is varying. We have shown 
that as long as the parameter e remains small enough, 
then there is a good adiabatic approximation describing 
the evolution of a Q-ball in an expanding background. In 
fact from the definition e = H/luq, we see that this de- 
creases with time in any universe with dominating fluid 
whose equation of state Wdom > — 1 (note this is not the 
same u as defining the Q-ball) . Hence in a matter and ra- 
diation dominated universe the adiabatic approximation 
improves with time. 



V. THE APPROXIMATION AND ITS LIMITS 

Using the change of variables we have introduced it is 
possible to analytically consider the behaviour of a Q-ball 
in a slowly expanding background. Let us analyse the 
approximations we have made in more detail. Consider 
the general potential given in equations 1)4. 9f) and (|3.7fl . 
This leads to an effective potential given by 



V,..ff = - | -<ur + — + - 



(5.11) 

An adiabatic evolution allows us to solve for Q-ball prop- 
erties using (|5.11|) and treating the scalefactor as a con- 
stant. Recalling the constraints arising in equation 12. 3|) 
for a static flat space solution, we can obtain instanta- 
neous constraints for the expanding Universe case: 
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We can also derive the scalefactor dependence of the Q- 
ball parameters. For example, using standard thin wall 
limit arguments j 1 1 1 ] we have the expression for the core 
field amplitude 
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and for the Q-ball radius 
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(5.13) 
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By introducing the dimension- 
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with the range < C < 1 we can write this as 
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(5.16) 



We see that these expressions, derived by considering the 
equations at a particular time, correctly predict the de- 
pendence of the parameters on the scalefactor. Similar 
calculations can be used to predict the background de- 
pendence of other, more complicated parameters. It is 
encouraging that the analytic estimates arising from the 
use of the adiabatic parameter e, and the thin wall ap- 
proximation lead to the correct values for the actual Q- 
ball parameters. Of course, as e increases, the validity of 
the analytic result comes into question, eventually break- 
ing down for large enough expansion rates. Fortunately 
though, it appears to hold well for a wide range of cosmo- 
logically interesting parameters. The breakdown of the 
approximation can be seen through the inhomogeneous 
charge distribution which emerges as e increases as seen 
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in figure 10) ■ These type of Q-balls are very interesting 
in themselves as they differ from the Q-balls discussed 
in with a homogeneous rotation frequency distribu- 
tion. In fact for significant expansion rates, they can 
evade the bounds imposed in 1(5.12(1 . which correspond 
to 12.3(1 . However, the bounds of ((5.12(1 do give us back- 
ground dependent constraints on the existence of nearly 
- homogeneous-w Q-balls. 

VI. CONCLUSIONS 

In this paper we have analysed Q-ball solutions in 
an expanding background. As expected, the presence 
of expansion changes the solution from the traditional 
form obtained in a Minkowski background, namely those 
which are spherically symmetric and have a constant ho- 
mogeneous frequency of rotation. However, by solving 
the system in conformal time, we can rewrite the equa- 
tions of motion for the Q-ball in an expanding back- 
ground in such a way that it mimics the case of the static 
background, the only difference being that the effective 
potential for the Q-ball given in equation 1(5.11(1 becomes 
explicitly time dependent as a result of the scale factor. 
This then allows us to introduce an adiabatic parameter 
e = with which to analyse solutions to the full equa- 
tions of motion. The smaller e is, the closer the evolution 
is to the original static spacetime case. Explicit calcula- 
tions of the deviations from static profiles, seen in Fig. 0] 
give us their dependence on e. This is important informa- 
tion when considering Q-balls in a cosmological context. 
Recalling that e is a function of both the expansion rate 
and the rotation frequency we see that in backgrounds 
where these effects play a role, even small variations in 
the frequency or the expansion rate can cause movement 



of charge. We numerically showed that Q-balls exist- 
ing in an expanding background conserved their charge, 
maintained a constant radius, whilst exhibiting oscilla- 
tions in their frequency, a reflection of the dependence 
on the expansion of the universe. The analytic approxi- 
mations we developed are able to recover the background 
dependence of these features as seen through equations 
((5.13|) and 1(5.16(1 . We are also able to place background 
dependent constraints on the existence of Q-balls with 
homogeneous rotation frequencies ((5. 12(1 . 

We have noticed in the simulations that the Q-ball so- 
lutions, whilst still existing in that they conserve their 
charge, now develop an inhomogeneous charge distribu- 
tion (see figure (0). Further study of the configurations 
could include predicting the general form of the charge 
density profile as a function of the background expan- 
sion. Given that the Q-ball charge profile is dependent 
on the expansion rate in a way which drives charge to 
the surface, we would anticipate that for large enough e, 
the charge is confined to the surface of the Q-ball. Such 
a configuration is unstable due to a large surface tension 
(defined by a positive energy-charge relation) and no bulk 
tension to hold it together. A similar situation for the 
case of non-topological cosmic strings is discussed in ^1 
and the same conclusions are reached. Any cosmological 
considerations of Q-balls where the oscillation frequency 
is of the order of the Hubble parameter should include 
discussions of stability against this kind of effect. This 
is particularly important in the early universe where we 
are dealing with rapid expansion rates. 
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